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EXHIBIT A 



^/ TECHNOLOGY CENTER R3700 

METHOD FOR ACQUIRING MRI DATA FROM A LARGE FIELD 
OF VIEW USING CONTINUOUS TABLE MOTION 

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH 

[0001] This invention was made with government support under Grant Nos. 
CA37993 and HL37310 awarded by the National Institute of Health. The United 
States Government has certain rights in this invention. 

CROSS-REFERENCE TO RELATED APPLICATIONS 

f00021This applicati on is based on U.S. Provisional Patent Application Serial 
No. 60/282.555 filed on April 9. 2001 and entitled "Method For Acnuirino MRI Data 
From A Large Field Of View Using Continuous Table Motion." 

BACKGROUND OF THE INVENTION 

f0QQ31fQOQ2] — The field of the invention is magnetic resonance angiography 
("MRA"), and particularly, studies of the human vasculature using contrast agents 
which enhance the NMR signals. 

ICL0.04] fOQ03] — Magnetic resonance angiography (MRA) uses the nuclear 
magnetic resonance (NMR) phenomenon to produce images of the human 
vasculature. When a substance such as human tissue is subjected to a uniform 
magnetic field (polarizing field SefJa), the individual magnetic moments of the spins in 
the tissue attempt to align with this polarizing field, but precess about it in random 
order at their characteristic Larmor frequency. If the substance, or tissue, is 
subjected to a magnetic field (excitation field Bi) which is in the x-y plane and which 
is near the Larmor frequency, the net aligned moment, MzM 7 . may be rotated, or 
"tipped", into the x-y plane to produce a net transverse magnetic moment M t . A 
signal is emitted by the excited spins, and after the excitation signal Bi is terminated, 
this signal may be received and processed to form an image. 

L0fl-0-51 fQQQ4J — When utilizing these signals to produce images, magnetic 
field gradients (G x G y and G z ) are employed. Typically, the region to be imaged is 
scanned by a sequence of measurement cycles in which these gradients vary 
according to the particular localization method being used. The resulting set of 
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received NMR signals, or-^views-^ are digitized and processed to reconstruct the 
image using one of many well known reconstruction techniques. 

I£2fifil{QG©§] — MR angiography (MRA) is the application of magnetic 
resonance imaging methods to the depiction of the human vasculature. To enhance 
the diagnostic capability of MRA a contrast agent such as gadolinium can be 
injected into the patient prior to the MRA scan. Excellent diagnostic images may be 
acquired using contrast-enhanced MRA if the data acquisition is properly timed with 
the bolus passage. 

f_Q_QLQZ]fQOO§] — The non-invasiveness of MRA makes it a valuable screening 
tool for cardiovascular diseases. Screening typically requires imaging vessels in a 
large volume. This is particularly true for diseases in the runoff vessels of the lower 
extremity. The field of view (FOV) in MR imaging is limited by the volume of the 
Sofia field homogeneity and the receiver coil size (typically, the FOV<48 cm on 
current commercial MR scanners). The anatomic region of interest in the lower 
extremity, for example, is about 100 cm and this requires several FOVs, or stations, 
for a complete study. This requires that the patient be repositioned inside the bore 
of the magnet, the patient be re-Iandmarked, scout images be acquired and a 
preparation scan be performed for each FOV. All of these additional steps take time 
and, therefore, are expensive. When contrast enhanced MRA is performed, the 
repositioning also necessitates additional contrast injections. 

£QQQ81f0007] — Recently gadolinium-enhanced bolus chase techniques have 
been reported which overcome this difficulty, K.Y. Ho, T. Leiner, M.H. de Hann, 
J.M.A. van Engleshoven, "Gadolinium optimized tracking technique: a new MRA 
technique for imaging the peripheral vascular tree from aorta to the foot using one 
bolus of gadolinium (abs). M Proc. 5th Meeting of ISMRM . p203, 1997. As described 
in U.S. Pat. Nos. 5,924,987 and 5,928,148, MRA data is acquired from a large field 
of view by automatically moving the patient table to a plurality of different locations 
during the scan and acquiring an image at each station. The movement of the table 
may be timed to follow the contrast bolus through the vasculature so that peak 
contrast is achieved at each station. 

IQflQgj fQQQfl — The result of prior bolus chase MRA methods is that one 
ends up with a plurality of images. These are manually or automatically registered 




with each other to provide a single image that covers the entire extended field of 
view. One difficulty with this approach, however, is that the separate images have 
different brightnesses and/or contrasts. As a result, there are discontinuities at the 
boundaries of images where they have been patched together. Another difficulty 
with the multi-station method is that valuable time is lost when the table is moved 
from one station to the next. During that time no image data is being acquired and 
further time is lost in bringing the spin magnetization into dynamic equilibrium before 
image data is acquired. In a three-dimensional scan this lost data acquisition time 
can mean that the scanning process does not keep pace with the moving contrast 
bolus and some image contrast is lost in later images. 

SUMMARY OF THE INVENTION 

[QQIQl tOOOQ] — The present invention is a method for acquiring an MR image 
from an extended field of view in which the patient table is continuously moved 
during the scan. The acquired MRI data is motion corrected to a common table 
reference position and a single, seamless MR image is reconstructed. MRI data 
acquired during the scan may be used to reconstruct images in real-time that enable 
the operator to monitor the progress of the contrast bolus and to adjust the speed of 
the table motion accordingly to insure peak contrast throughout the scan. 

fQQj_lTf0046] — A general object of the invention is to produce a single MR 
image over a field of view which exceeds the normal field of view of the MRI system. 
The subject is moved continuously through the MRI system and views are acquired 
along with table position data. Each view is position corrected using the associated 
table position data to produce a single array of MRI data which is used to 
reconstruct an image. 

r0dl2] fQQ44^ — Another object of the invention is to better time the 
acquisition of MRI data with the movement of the contrast bolus through the subject. 
Once the scan is begun no time is wasted. For a given image resolution and pulse 
sequence, the table is moved continuously at the maximum possible velocity. Table 
velocity can be slowed as needed to remain properly aligned with peak contrast as it 
traverses through the subject. 



IQQ±3\[QQ4Z] — The foregoing and other objects and advantages of the 
invention will appear from the following description. In the description, reference is 
made to the accompanying drawings which form a part hereof, and in which there is 
shown by way of illustration a preferred embodiment of the invention. Such 
embodiment does not necessarily represent the full scope of the invention, however, 
and reference is made therefore to the claims herein for interpreting the scope of the 
invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 

10014} fQ&43] — Fig. 1 is a block diagram of an MRI system which employs 
the present invention; 

JQQ±5]f0044^ — Fig. 2 is a graphic representation of a pulse sequence 
performed by the MRI system of Fig. 1 to practice a preferred embodiment of the 
invention; 

[PP16] fQQ4£] — Fig. 3 is a schematic representation of an exemplary imaging 
volume within the bore of the MRI system of Fig. 1 ; 

[QJQ1X[{Q&4§] — Fig. 4 is a schematic representation of a data matrix into 
which data acquired with the pulse sequence of Fig. 2 is stored; and 

JflPiai [0047] — Fig. 5 is a flow chart illustrating the steps employed in 
practicing the preferred embodiment of the invention. 

GENERAL DESCRIPTION OF THE INVENTION 

[0Q_19j {OG4&j — The present invention allows the subject to be imaged over a 
field of view (FOV) which can be many times larger than the static FOV allowed by a 
typical MRI scanner. The invention allows the MRI scanner to run uninterrupted 
while the patient table is continuously moved through the portion of the scanner 
volume which presents the best combination of homogeneous static magnetic field 
and linear magnetic field gradients. This region is typically known as the scanners 
-Ilsweet spot.-I The invention allows the acquisition of arbitrarily large FOVs along 
the table motion direction with one uninterrupted data set. Two and three- 
dimensional imaging are possible. The method of imaging during continuous table 
motion can be coupled with other techniques, including real-time MR imaging, to 
provide real-time control of the table motion and scan parameters. The invention is 



primarily directed toward MR angiography, but is also a general MRI imaging 
method. 

[Qfl2fl] fQQ49] — The invention can be described by referring to Fig. 3 which 
schematically shows the bore 10 of an MRI scanner. For this discussion a slab 12 
with coronal orientation is being imaged. The readout direction (X axis) is along the 
direction of motion of the table and the slab select direction is along the Z axis. For 
2DFT imaging, phase encoding is performed along the Y axis and for 3DFT imaging, 
phase encoding is performed along both the Y and Z axes. 

\QQ2i] \QQ2&\ — We define FOVtot as the full desired field of view along the X 
direction and POVxEQ^x as the field of view along X for the readout of any single 
NMR signal. We assume that FOVtot > FOVxEQV* thus necessitating the invention 
described here to provide an image for the full FOV if continuous table motion is 
desired . 

\QQ22l\QQ2A] — As the subject is moved through the scanner along the x-axis, 
MRI data is acquired from the slab 12 by a series of imaging pulse sequences. 
During each pulse sequence the readout of the encoded k-space data is performed 
in the X direction. The raw k-space data is typically anti-alias filtered and 
subsequently sampled at a bandwidth assuming some PQV yFQVv in the readout (X) 
direction. 

\QQ231\QQ2Z) — Referring to Fig. 4, the raw k-space data are acquired in a 
series of imaging pulse sequences, a Fourier transformation (FT) is performed on 
each NMR signal in the readout direction and the result is placed in a matrix of 
memory locations 16 which represents the desired resolution in the X, Y and Z 
directions. This once-transformed view is a line of data placed in the matrix 16 and 
shifted in the X direction to a position that represents the patient position at the time 
the NMR echo signal was acquired. The position of the patient table relative to a 
reference location is known for every echo readout. The data acquisition is 
performed until all of the data matrix 16 is filled and the table has reached the final 
position. The final image is formed by Fourier transforming along the remaining Y 
and if 3D, the Z directions. 

fflQ241 IhfiLMtow 

atMtraj^QQsitioninq of thajaliepLan^ acquired NMR echo 
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signal. The MR signal measured durincia_s_taa^^ 
cLes_ccib_e_dJ)y ; 



[0023] 



Th e fo ll owing i s a d i scuss i on of th e th e ory of signa l a l t e ration du e to 



arb i trary pos i t i on i ng of tho pati e nt and th e adjustm e nt of the acqu i red NMR ocho 
s i gna l . Tho MR s i gna l measured dur i ng a standard i mag e acqu i s i tion can bo 
d e scr i b e d by: 



where x is the frequency-encoding (readout) direction, y is the phase-encoding 
direction, GxQ* and G n are the respective gradients, m(x,y) is the transverse 
magnetization of the object in spatial coordinates, ^ is the gyromagnetic ratio, and t y 
is the duration of the y gradient. The index n is the repetition number of the 
acquisition and runs from 0 to N-1 where N is the total number of phase encodings 
along the Y axis. 

[002 4 ] — An obj e ct, shift e d from its in i t i a l r e f e r e nc e pos i tion along x som e 
arb i trary distanco has pos i tion An object, shifted from its initial reference position 
afongj^Sjpme^ Tho s i gnal from th i s sh i fted objoct is 

th e n giv e n by: 



$1,(0 = J \m(x>y)e 




dxdy 



(D 



Now by subst i tution of var i ables 



w e hav e 

y 



S„(t) = jjm(x - A,y)e~ irG " yly e~ i>G * x 'dxdy 

Ro-W-by_s.u bstitution oJLv_apables 
x' = x - A and dx' = dx 
we have 

S ' n (t) = \\m{x< ^e^'e^^^dx* dy 

crrrcr 



and 



y e ~ iyG * xU 



dx* dy 



where ^ is a phase term representing the positional shift of the object in k- 

space. It is assumed that Gn and A are known for any specific phase encoding. 

+iyG x bt 

The phase factor required to recover or unshift the signal Sn(t) is then e 

[0025] \QQ2S] — There are two possible methods to position each echo using 
the previous theory. The direct method in which the shifted data is repositioned by 
applying the phase factor to the continuous presampled, raw k-space data, and the 
time and memory efficient method in which a combination of phase shifting and 
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positional shifting is applied to the sampled data. These are described in more 
detail in the description of the preferred embodiment of the invention. 

DESCRIPTION OF THE PREFERRED EMBODIMENT 

tOJQ26J fOQ2^ — Referring first to Fig. 1 , there is shown the major components 
of a preferred MRI system which incorporates the present invention. The operation 
of the system is controlled from an operator console 100 which includes a keyboard 
and control panel 102 and a display 104. The console 100 communicates through a 
link 116 with a separate computer system 107 that enables an operator to control 
the production and display of images on the screen 104. The computer system 107 
includes a number of modules which communicate with each other through a 
backplane. These include an image processor module 106, a CPU module 108 and 
a memory module 113, known in the art as a frame buffer for storing image data 
arrays. The computer system 107 is linked to a disk storage 111 and a tape drive 
1 12 for storage of image data and programs, and it communicates with a separate 
system control 122 through a high speed serial link 115. 

JHQ271 [0027] — The system control 122 includes a set of modules connected 
together by a backplane. These include a CPU module 119 and a pulse generator 
module 121 which connects to the operator console 100 through a serial link 125. It 
is through this link 125 that the system control 122 receives commands from the 
operator which indicate the scan sequence that is to be performed. The pulse 
generator module 121 operates the system components to carry out the desired 
scan sequence. It produces data which indicates the timing, strength and shape of 
the RF pulses which are to be produced, and the timing of and length of the data 
acquisition window. The pulse generator module 121 connects to a set of gradient 
amplifiers 127, to indicate the timing and shape of the gradient pulses to be 
produced during the scan. The pulse generator module 121 also receives patient 
data from a physiological acquisition controller 129 that receives signals from a 
number of different sensors connected to the patient, such as ECG signals from 
electrodes or respiratory signals from a bellows. And finally, the pulse generator 
module 121 connects to a scan room interface circuit 133 which receives signals 
from various sensors associated with the condition of the patient and the magnet 



system. It is also through the scan room interface circuit 133 that a patient 
positioning system 134 receives commands from the pulse generator module 121 to 
move the patient through the scanner to perform the scan in accordance with the 
present invention. The current position of the table at any time during the scan is 
read into the system control 122 and is used to adjust the acquired NMR data 
according to the present invention as will be described in more detail below. The 
operator can control the operation of the patient positioning system 134 through the 
keyboard and control panel 102. This includes controlling the velocity of table 
motion during the scan. 

[0.028] [0028] — The gradient waveforms produced by the pulse generator 
module 121 are applied to a gradient amplifier system 127 comprised of G^, G y 
and GZQz amplifiers. Each gradient amplifier excites a corresponding gradient coil 
in an assembly generally designated 139 to produce the magnetic field gradients 
used for position encoding acquired signals. The gradient coil assembly 139 forms 
part of a magnet assembly 141 which includes a polarizing magnet 140 and a whole- 
body RF coil 152. A transceiver module 150 in the system control 122 produces 
pulses which are amplified by an RF amplifier 151 and coupled to the RF coil 152 by 
a transmit/receive switch 154. The resulting signals radiated by the excited nuclei in 
the patient may be sensed by the same RF coil 1 52 and coupled through the 
transmit/receive switch 154 to a preamplifier 153. The amplified NMR signals are 
demodulated, filtered, and digitized in the receiver section of the transceiver 150. 

[0029] — Tho transm i t/rece i v e s w i tch 154 i s contro lle d by a 6i gna l from 
tho pu l so g e nerator modulo 121 to o l octr i ca l ly connocttho RF amp li f i er 151 to tho 
coi l 152 dur i ng th e transmit mod e and to conn e ct th e pr e amp li f ie r 153 dur i ng tho 
r e c ei v e mod e . Th e transm i t/r e c ei v e sw i tch 154 a ls o e nab le s a s e parat e RF loca l 
co il to bo usod dur i ng tho rocoivo mod e . I n th e pr e f e rred e mbodim e nt of tho 
i nv e nt i on a s i ngl e , stat i onary l oca l co il doocr i bod be l ow and shown in Fig. 7 i s 
sw i tch e d i nto op e rat i on. I n th e a l t e rnat i v e , a mu l t i- s e gm e nt l oca l co i l doscr i bod 
b el ow and shown i n Fig. 8 may also b e us e d. I n th i s cas e , th e s e parat e s e gm e nts of 
tho l ocal co i l are sw i tchod i nto op e ration by th e transm i t/r e co i vo sw i tch 15 4 as th e 
pat ie nt i s trans l at e d to succossivo stat i ons during tho scan. 



[0029]The NMR signals picked up by the RF local coil are digitized by the 
transceiver module 150 and transferred to a memory module 160 in the system 
control 122. Whon tho scan i s compl e t e d and an e nt i r e array of data has b ee n 
acqu i red in tho momorv modu l o 160. an An array processor 161 operates to Fourier 
transform the data into an array of image data. This image data is conveyed 
through the serial link 1 15 to the computer system 107 where it is stored in the disk 
memory 111. In response to commands received from the operator console 1 00, 
this image data may be archived on the tape drive 1 12, or it may be further 
processed by the image processor 106 and conveyed to the operator console 100 
and presented on the display 104. 

[0030] For a more detailed description of the transceiver 150, reference is 
made to U.S. patent Nos. 4,952,877 and 4,922,736 which are incorporated herein by 
reference. 

[0031] While many pulse sequences may be used to practice the present 
invention, in the preferred embodiment a 3D gradient-recalled echo pulse sequence 
is used to acquire the NMR data. Referring particularly to Fig. 2, an RF excitation 
pulse 220 having a flip angle of 45-f, is produced in the presence of a slab select 
gradient pulse 222 to produce transverse magnetization in the 3D volume of interest 
as taught in U.S. patent No. 4,431 ,968. This is followed by a phase encoding 
gradient pulse 224 directed along the z axis and a phase encoding gradient pulse 
226 directed along the y axis. A readout gradient pulse 228 directed along the x axis 
follows and a partial echo (60%) NMR signal 230 is acquired and digitized as 
described above. After the acquisition, rewinder gradient pulses 232 and 234 are 
applied to rephase the magnetization before the pulse sequence is repeated as 
taught in U.S. patent No. 4,665,365. 

[0032] As is well known in the art, the pulse sequence is repeated and the 
phase encoding pulses 224 and 226 are stepped through a series of values to 
sample the 3D k-space in the field of view. As w ill b e descr i b e d i n mor e deta il 
b el ow, i t i s th e ord e r i n which th e s e phas e e ncod i ng grad ie nts ar e st e pped, and 
h e nc e th e ord e r i n wh i ch k spac e i s sampl e d, wh i ch i s th e subj e ct of th e pr e s e nt 
i nv e ntion. In the preferred embodiment 32 phase encodings are employed along the 
z axis and 256 phase encodings are employed along the y axis. Sampling along the 



kx axis is performed by sampling the echo signal 230 in the presence of the readout 
gradient pulse 228 during each pulse sequence. It will be understood by those 
skilled in the art that only a partial sampling along the k x axis i smav be performed* 
and if this is done, the missing data is computed using a homodyne reconstruction 
or by zero filling. This enables the echo time (TE) of the pulse sequence to be 
shortened to 1 .0 ms and the pulse repetition rate (TR) to be shortened to 4.5. 4 .5 ms. 

[0033]The preferred embodiment of the invention is a 3D CE-MRA scan of a 
subject after injection of a contrast agent. When the scan prescription is entered by 
the operator, scan parameters such as the imaging pulse sequence FOV and 
resolution are defined along all axes. Unique to the present invention, the FOV to t is 
also defined. Based on this information a data matrix 16 is defined as illustrated in 
Fig. 4. Referring particularly to Fig. 5, after the injection of the contrast agent at 
process block 240, 2D images are rapidly acquired and reconstructed to display the 
vasculature at the starting boundary of the FOVtot. Arrival of the contrast bolus is 
detected at process block 242 by observing when the arteries brighten. At this time 
table movement is initiated as indicated at process block 244. There are other 
methods for detecting bolus arrival, including automated methods such as that 
described in U.S. Pat. No. 6,167,293. 

[0034] A loop is then entered at 246 in which MRI data is acquired as the 
table moves the patient through the sweet spot of the scanner. A table reference 
location is also established at this time. As indicated at process block 248, a view is 
acquired by performing the above-described pulse sequence. Many different view 
orders may be employed and the prescribed view order is stored in table 250. After 
acquisition of the NMR echo signal, the current table position is also acquired, as 
indicated at process block 252, and the NMR echo signal is Fourier transformed. 
along the x-axis as indicated at process block 254. As indicated at process block 
256, the proper location in the data matrix 16 for the acquired view is then calculated 
and the data is stored therein as indicated at process block 258. 

[0035]The location for each acquired view in the data matrix 16 is determined 
by two factors, the view number and the location of the patient table at the time the 
view was acquired. This is illustrated in Figs. 4 and 6 for a two ^dimensional 
acquisition with a monotonic view order. The first view is acquired at echo location 



x(0) which serves as a reference location. Subsequent views are acquired as the y 
axis phase encoding is stepped through its 128 values. The location in the data 
matrix 16 along its k y axis is determined by the G *G y phase encoding for the view. 
As indicated by arrow 4SQ t260. the center of each successively acquired NMR echo 
signal after Fourier transformation along x is also shifted along the x-axis due to 
table movement. The amount of this shift from the reference position x(0) is 
measured at the time of view acquisition using an electronic spatial encoder or by 
using the following equation: 



whoro Vr e f i s tho tab l o ve l oc i ty, t is th e tim e past sinc e th e b e ginn i ng of tho scan, n 
i s th e numb e r of pu l s e s e qu e nc e r e p e t i t i ons sinco tho start of tho scan and TR i s th e 
t i m e r e qu i r e d for on e pu l s e s e qu e nce r e p e t i t i on. 

x(n) = V ref -t 



x(n) = V ref -n-TR 

[0037] — where Vmf isJtheJtahl^^ of the 

scan, n is thejoumfa^^ 

TR is the time requi rejlfcjLia^ repetition. The storage pattern is 

similar to a usual Cartesian k-space sampling pattern, but it is skewed by the table 
motion. When all the phase encoding views have been acquired, the process 
repeats and each new view of once transformed data is concatenated to the 



or a l t e rnat i v el y 
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previous acquired kv]& view. It should be apparent that in a 3D acquisition a series 
of kz)& phase encoded views are al§Q-acquired at each k y phase encoding, and 
these are shifted and stored in a 3D data matrix 16 in the same manner. 

[Qfl3Sl fQQ3gl I t can bo apprec i ated tha t ln general, the table motion might 

be such that the displacement of an acqu i r e d v ie w along tho x ax i s^ from one 
repetition to the next mayis not be-an integral number of ^pixels. The data for such 
a-vie wthese repetitions must be shifted with a precision of less than one pixels 
Th e r e ar e two possibl e m e thods w i th e qu i val e nt outcom e s that can b e app lie d to 
perform this data shift. Tho f i rst ca ll ed th e d i r e ct m e thod app l i e s on e phase shift to 
tho raw, k spac e data to sh i ft i t b v to maintain consistency. The displacement X(n) 
is known and can be split into two parts as follows: 

X(n^=m - A . Th e s e cond method ca ll ed tho t i mo and memory e ff i ci e nt 
m e thod uses a phys i ca l sh i ft of th e onc e Four ie r transform e d data p l us an add i t i ona l 
phas e man i pu l ation of tho sampl e d raw k spaco data. 

rOQ-3Z U 0039] — At Tho d i rect mothod . S » (t) is rocovor e d from S' ^ftj-by 
apply i ng th e phas e t e rm during d e modu l at i on of th e raw cont i nuous NMR signa l . 
Tho process of demodu l ation of th e raw signa l as w ell as th e app li cat i on of phas e 
accrua l dur i ng th e d e modulation is wo ll known to those studi e d in th e art. Th e 
corr e ct e d NMR data i s th e n samp le d to giv e th e samo number of data po i nts as aro 
r e qu i r e d for th e numb e r of p i x el s r e pr e s e nt i ng the FO\Aai i n th e X dir e ction x + 
^(e xamp le Ntoi ~ 2CM8). After Fourior transform i n X, th e onc e transform e d hybr i d 
spac e signal i s H(x, k y ). Tho r e g i on roprosonting FOV A i s reta i n e d by tak i ng th e 
product of H(x,k y ) and a window function W(x A) wh i ch has b ee n sh i ft e d by a 
d i stance n ) (5) 

^fa&L^A. Th e simp le st w i ndow funct i on i s th e r e ctang l e funct i on Roct(x), but othor 
w i ndow funct i on s wh i ch ar e non zero outs i d e FOVx ar e possib l o. I t is assumed that 
v i a th e convo l ut i on th e or e m that convo l ut i on of th e Fouri e r transform of W(x) w i th 
Sn (t) is equ i va l ent to tho product of W(x) and H(x,k y ) aft e r Four ie r transform of S n {t)r~ 
Succ ee d i ng signals are added to the hybr i d soaco matr i x x is the pixel to pixel 



-JAz 



Ax = 



FOV, FOV, r 



N. 



N„ 



where Nv and N,n> are the number of pixels in the X direction for the echo readout 
apd_fo r th e ir sp e o i fio kv va l u e s gftor w i ndow i ng to th e ir sp e oif i o known the total FOV 
respectively. Also, m is the largest integer such that m-A and FOV^x<Xfnl and F>(n\ 
is the sub-pixel-sized remainder which forces equation (5^ to hold. As hefore. n is 
the index on the pulse sequence repetition number. The actual displacement 
positioning consists of two steps, coarse and fine, presented here in reverse order of 
application for_ease_of underslaridjng . 

[0040] — St T i m e and memory e ffic ie nt m e thod . Th e signa l S' B {t)-is- 
samp le d for a spoo i f i o FOV* at somo arb i trary FOW. Th e number of samples (N*) i o 
le ss than th e numb e r of p i x el s r e present i ng th e f i nal imag e (Ntot) w i th FOVt e t i n th e X 
direction. As an oxamplo Nt e t ~ 2048, N» ~ 256. Each pixo l w i dth roprooonto a 
d i stanc e (Dp i *) on th e obj e ct wh e r e 



Actua l posit i on i ng of data i n th e Hybrid matrix roqu i res two stops. Thoy aro 
pr e s e nt e d i n r e v e rs e ord e r of app li cation for e as e of und e rstand i ng. 

IM381{0©44} — Step 2, coarse positioning: Aft e r th e sampl e d k spaco data 
\&Ih& Fourier transformed a l ong x i t can bosional is shifted hv m pixels as given in 
the first term of equa tion (5) and directly placed miflto, the hybrid matrix 16 with a 
positional accuracy of *Bpi^. The data then requires additional sub-pixel correction, 
bv the amo unt S(n) to maintain data consistency. This process is performed in step 
1. 

IQQ29H QQ42\ — Step 1 , fine positioning: Prior to Fourie r transform. 
transformation in X, a phase tefmiwjsl representing the sub-pixel correction 5m ^ is 
applied toalono the sampled data. Th e subpixo l sh i ft (A') is g i v e n by: 



A' - A n > D ftm 



whoro n io chosen such that A' < D^. Th e samp le d s i gna l w i th this fin e corr e ct i on h 
th e n g i von by: 



scho: i.e. 



_S n (k x ,k y (n)) = e N * .^(k^k^n))^ ^ 

Fourier transformation of S n along the number of tho samp l ed data po i nt rang i ng 
from 0 to Nx - 1 . Th e sign app lie d to th e phas e t e rm is d e p e nd e nt on whothor tho 
coars e corr e ct i on w ill b e a part i a l pix el l e ss than or gr e at e r than th e actua l pos i t i on 

Ax-axis followed bv the coarse positioning of Step 2 yields the desired S " ( x > k y >( n ))_ 
In principle it would be possible to perform the entire d isplacement positioning using 
the Phase manipulation approach of eouation (7Y but because multiple Fourier 
transformations of length Ntnt would be required for each kv value a cquired, this 
"direct method" is not preferred because the two part position shifting and phase 
shifting has been found to be computationally more efficient. 

£004Q]{004S] — It can be seen that after one complete cycle of k y phase 
encoding values has been completed and stored in the data matrix 16, the process 
repeats itself and the phase encodings are applied a second time. Depending upon 
the number of phase encodings and the patient table velocity, the second time a 
view is acquired for a specific phase encoding, the placement of the echo signal 
samples in the data matrix 16 along the x direction may or may not overlap with 
some of the echo signal samples acquired for the previous measurement at that 
phase encoding. In the 2D example shown in Fig. 4, such a view is shown at x(n) 
and the timing is such that there are no overlapping samples with the previously 



acquired view at x(0). In general some overlap will occur and the redundant, 
overlapping signal samples are averaged to improve SNR. 

fQ_Q4_1_]fG044] — Referring again to Fig. 4, views are acquired and stored in 
the data matrix 16 until samples have been acquired for the entire extended field of 
view FOVtot. The last view is indicated at 262 in Fig. 4, and after it is acquired and 
stored as determined at decision block 264, the scan is complete and a single image 
is reconstructed at process block 266 using the samples in data matrix 16. In the 3D 
acquisition of the preferred embodiment, this reconstruction includes Fourier 
transformation of the data matrix 16 along the k y and kz axes and calculation of pixel 
intensities from the resulting complex numbers. Typically, a two-dimensional 
projection image is then produced from the resulting three-dimensional array of 
image pixel intensities. The well-known maximum intensity pixel (-^MIP-^) technique 
is preferred for producing the projection image and a pre-contrast mask image may 
also be subtracted to suppress signals from non-vascular tissues. 

10042] fQQ4£| — As shown in Fig. 4, the extended field of view FOVtot over 
which complete k-space sampling is acquired is smaller than the x axis extent over 
which samples are acquired. The regions 268 and 270 of partial sampling at the 
beginning and end of the scan can also be reconstructed to extend the total field of 
view, but it can be appreciated that image quality will decline as fewer samples are 
used in the reconstruction. 

IQM31{O046] — Referring again to Fig. 5, after each view is acquired and 
stored during the scan, data in the data matrix 16 may be employed to reconstruct a 
two-dimensional monitor image as indicated at process block 272. Real-time 
images are thus produced for the operator who can use the information therein to 
control the scanning process. In addition to changing the usual scan parameters, 
the operator can adjust the table velocity as indicated at process block 274. This 
might be done, for example, to better match the speed at which the contrast bolus is 
moving through the extended field of view FOVtot. 

[0044} [OCM 7] — The real-time 2D monitor image may be produced in a 
number of ways. The preferred method is to reconstruct a full maximum intensity 
pixel (MIP) projection through the image stack as the data is being acquired. Unlike 
typical static 3D image acquisitions, a partial 3D data set can be reconstructed as 



-JZz 

the hybrid matrix fills. For Nz slices and Ny phase encodes in Y, Ny Fourier 
transforms in the Z direction and Nz Fourier transforms in the Y direction can be 
performed with each TR. Each set of these described Fourier transform sets in Z 
and Y (including the previous transform in X) provide a partial 3D image set that is of 
one pixel length in the X direction. Ny pixels in the Y direction and Nz pixels in the Z 
direction. The maximum pixel value of a projection through this partial image set as 
well as earlier partial data sets can be found at the same time. The projection may 
be in the Z direction directly or in oblique angles through previously reconstructed 
partial image sets. The 2D MIP images can be displayed on the real-time monitor 
as each portion of the 3D data set is filled. Also, specific individual slices from the 
3D data set can be viewed in real time rather than the MIP projection. 

£QQ451fQQ4S] — Alternatively the same process can be performed before the 
hybrid matrix is completely filled. No matter what phase encoding strategy is used, 
the phase encode order will periodically return to the center of k-space during the 
acquisition. Images viewed would possibly have a lower spatial resolution (in Y and 
Z) but would be more temporally recent than allowing full hybrid space filling. Th e s e 
r e a l -tim e 2D i mages may bo usod to mon i tor scan progr e ss. Th i s not on l y i ncludes 
th e physica l pos i t i on of th e acquis i t i on on tho pat ie nt at that t i mo but can b e usod to 
monitor th e contrast bo l us i ts e lf. Th e bolus t i mo pos i tion i nformat i on can b e usod to 
r e duc e or poss i b l y i ncr e as e th e ve l oc i ty of th e scan dur i ng tho runoff study. 

r.0.0461 These real-time 2D or low resolution 3D images may be used to 
jODQOJlQLJScaDj aEDC ^ acq uisition 
pjxihejaMei^ 

toiiisj^ information can JaejJsedJcLt ^^ 

Yelojsi^ 

[0047]{0©4©] — The table velocity to be used for the continuous motion is an 
independently chosen parameter. There are guidelines for its selection. In order to 
match table motion to the playout of phase encoding views the table should move a 
complete field-of-view (EOV^EQ^) in the time required for application of one 
complete cycle of the phase encodings. At this velocity, consecutive measurements 
of the same view will abut and have neither an overlap ornor an intervening gap 
along the readout x direction. We define this as the -preference table velocity" and 



for 2D imaging is given by: 



v FOV x cm 

refaD- N R se( . 



where N y is the number of y phase encodes and TR is the repetition time of the 
imaging pulse sequence. For 3D imaging in order to allow for both phase encode 
directions the reference table velocity is given by: 



FOV r 



ref * D N y .N 2 .TR 



cm 



sec 



where is the number of z phase encodes. Note there is no dependence on the 
chosen resolution in X, only the chosen readout EQV xFOVv . The table can be 
slowed or even stopped if redundant data are desired. In that case the redundant 
data can replace or be averaged with the non-zero data already stored in the data 
matrix 16. 

J£Q482{GQ£0] FQV>£Q^ can be variable during the scan. An acquisition 

can be performed with higher table velocity and thus greater temporal efficiency if a 
portion, or all of the hybrid space is acquired with a larger FOV x . EQV *FOVv can 
also be reduced if a higher spatial resolution is desired for a limited region within the 
FOVtot. Regardless of changes in table speed or readout FOV x , a single, data matrix 
16 is formed from which a seamless image of the extended field of view FOVtot can 
be reconstructed. 

fQQ491We have shown the feasibility of acquiring 2D and 3D images whir.h 
cover longitudinal fields-of-view many times lo nger than the readout FOV possible in 
a single static image accuisition. The data for these images are acquired in 
conjunction with continuous tahle motion. This approach eliminates discontinuities 
in imaoe signal at station boundaries as well as the ne ed for overlap in FOV The 
total scan times reouired for the moving table MR A acquisition vs. a ge neric fixed 
multi-station approach can be compared. For equal readout FOVs the continuous 



motion approach of the present invention provides reduced overall scan time 
because time dedicated to station-to-station motion and dum my repetitions is 
eliminated. 

r00501 With the introduction of shorter bore, less claustrop hobic MRI systems. 
the present invention provides the ability to image large longitudinal fields of view in 
such systems. In addition, hv using one of the new fast 3D steady state pulse 
seouences the invention allows whole body screening with one fast exam. This can 
be done using targeted contrast aoents which determine, for examnle the extent of 
metastas es for certain cancers or the quantitative assessment of systemic thromus. 

r00511ln summary, continuously moving table acouisitions have been 
demonstrated for 7DFT and 3DFT MR acquisition using gradient e cho and spin- 
echo methods. Phantom and in vivo im ages are sharp and appear devoid of motion 
artifacts as well as interference effects between consecutive FOVs. The method 
offers the potential for high resolution 2D and 3D MRI of extended objects. 
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CLAIMS 



1. In a magnetic resonance imaging (MRI) system having a defined field 
of view (FOV), a method for producing an image of a subject over an extended field 
of view (FOVtot) which is larger than the FOV. the steps comprising: 

a) moving the subject through the MRI system such that the extended 
field of view (FOVtot) passes through the defined field of view (FOV); 

b) continuously acquire NMR data from the subject as it is moved through 
the FOV by repeatedly performing an imaging pulse sequence which acquires NMR 
data comprising a view of the subject; 

c) adjusting each view acquired in step b) using subject position 
information; 



d) storing each adjusted view in a data matrix; and 

e) reconstructing an image using the data matrix. 



2. The method as recited in claim 1 in which the MRI system has a table, 
and step a) is performed by: 

i )plac i ng th e subj e ct on tho tablo; and 

ID moving the table. 



3. The method as recited in claim 2 in which the table is moved 
continuously while performing step a). 



5. The method as recited in claim 2 in which the table is moved at 
different velocities while performing step a). 



6. 
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2± The method as recited in claim 2 which includes: 
injecting the subject with a contrast agent: and 

§r Tho mothod as rooit e d i n c l aim 2 wh i ch i nc l udes: i nj e cting th e subject 
w i th a contrast ag e nt; and 

in which the table is moved at a velocity which tracks the contrast agent as it 
moves through the extended field of view (FOVtot). 

L The meth ojj_as_Jiecjte d in claim 5 which includes: 

& Th e m e thod as r e cit e d i n c l a i m 5 wh i ch i nc l ud e s: reconstructing 

monitoring images from th e data matr i x during the performance of step a) from data 

stored in the data matrix. 

9. The method as recited in claim 1 in which step c) includes adjusting the 
location in the data matrix in which the view is stored in step d). 

11. The method as recited in claim 1 in which step c) includes adjusting the 
phase of the NMR data in the view. 

13. The method as recited in claim 2 in which step c) includes adjusting the 
location in the data matrix in which the view is stored in step d) as a function of the 
table location at the time the view is acquired in step b). 

15. The method as recited in claim 2 in which step c) includes adjusting the 
phase of the NMR data in the view as a function of the table location at the time the 
view is acquired in step b). 



The methoc Las_xe_cite_d_ in claim 2 in v\riiLC.h_sjep^^ 

D acquiring tabJteJLocation information as each view is acquire( 

and 
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44^ Th e m e thod as r e c i t e d i n c l aim 2 i n wh i ch st e p b) further i nclud e s: 

i ) acqu i r i ng tab le l ocat i on i nformation as e ach v i ew i s acqu i r e d; 

the table location information is used in step c) to adjust each corresponding 

view. 

JJL The method as recited in cl aim 1 1 in which st ep c) includes: 

D P-ej±ojjningja_Ea data in the view: 

and 

43t Th e m e thod as r e c i t e d in c l a i m 11 i n wh i ch st e p o) inc l udes: 

t) p e rform i ng a Four ie r transformation of th e NMR data i n th e view; and 

ii) calculating a location in the data matrix for the transformed view as a 
function of the table location at the time the view was acquired in step b). 

JJL The method as recited in claim 1 in which st ep c\ includes: 
0 adjusting the phase of the NMR data in the view: 
H) Fourier transforming the phase adjusted NMR data in the view: and 
iH} adjusting the locatio n in the data matrix in which the Fourier 
transformed view is stored in step d^ as a function of subject location at the time the 
view is acquired in step b^ with respect to a subject refe rence location 

2L The method as recited in claim 1 in which the performance of the 
imaging pulse sequence in step b^ includes: 

i) producing a readout magnetic field gradient during the acnuisition of 
said NMR data comprising a view, and the readout mag netic field gradient is 
oriented in the same direction as subject movement. 

ZL The method as recited in claim 1 4 in which step c) includes: 

1} Fourier transforming the acquired view: and 

ii} adjusting the location in the data matrix in which the Fourier 
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transformed view is stored in step d^ as a function of suhjent location at the time the 
view is acquired in step with respect to a subject reference location. 



2± The method as recited in claim 1 in which the data matrix 
dimensional array of data. 



22= The method as recited in claim 1 in which the data matrix is a three- 
dimensio nal array of data 

28. 



2SL In a magnetic resonance imag ing (MRW s ystem, the improvement 
comprising: 

a} a table for supporting a subject and for moving the subject through a 
defined field of view (FOV) of the MRI system: 

a pulse generator for operating the MRI system under the direction of a 
pulse sequence to continuousl y acguire a series of NMR data views of the subject 
as the subject is moved through the FOV: 

fi) means for adjusting each acquired view as a function of subject 
location at the time the view is acquired with resnect to a reference suhiect location: 

dj a memory for storing the adjust ed views as a data matrix: and 

a) means for reconstructing an image from data in the data matrix which 
has a field of view in the direction of table motion which is larger than the defined 
FOV. 

m 



The improvement as recited in claim 18 in which element r.\ includes: 
D means for Fourier tra nsforming each acquired view: and 
li) means for storing the Fourier transformed view in the data matrix at a 
location determined bv the subject locatio n at the time the view was acquired 

23* The improvement as recited in claim 18 which also includes: 
f) means for reconstructing an imaoe from data in the data matrix as the 
subject is moved through the defined FOV and views are heino acquired 

1L The improvement as recited in claim 20 which also includes: 
gj means for controlling the velocity of table motion as views are being 
acquired. 



2£L 
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ABSTRACT OF THE DISCLOSURE 

MRA data is acquired from a large region of interest by translating the patient 
through the bore of the MRI system as a three-dimensional MRA data set are 
acquired. Patient table movement is controlled to track a bolus of contrast agent as 
it passes through the region of interest. Fluoroscopic images may be acquired 
during the scan to enable accurate bolus tracking. A seamless image of the entire 
region of interest is reconstructed. 



